Cold streams in early massive hot haloes 
as the main mode of galaxy formation 
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c^. The massive galaxies in the young Universe, ten bilhon years ago, formed stars 

at surprising intensities.^'^ Although this is commonly attributed to violent 
^ ! mergers, the properties of many of these galaxies are incompatible with such 

rn ! events, showing gas-rich, clumpy, extended rotating disks not dominated by 

I spheroids. ^'■^'^ Cosmological simulations^ and clustering theory^ ' are used to 

O ' explore how these galaxies acquired their gas. Here we report that they are 

00 ' stream-fed galaxies, formed from steady, narrow, cold gas streams that pene- 

^ ■ trate the shock-heated media of massive dark matter haloes.®'^ A comparison 

O ■ with the observed abundance of star-forming galaxies implies that most of the 

^ ! input gas must rapidly convert to stars. One-third of the stream mass is in gas 

^ I clumps leading to mergers of mass ratio greater than 1:10, and the rest is in 

I smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies 

forming stars at a given rate are fed by smooth streams. The rarer, submillime- 
tre galaxies that form stars even more intensely^'^^'^^ are largely merger-induced 
starbursts. Unlike destructive mergers, the streams are likely to keep the ro- 
tating disk configuration intact, although turbulent and broken into giant star- 
forming clumps that merge into a central spheroid. ^^'^^ '^'^^ This stream-driven 
scenario for the formation of disks and spheroids is an alternative to the merger 
picture. 



Star-Formation Rate versus Halo Growth Rate 

It appears that the most effective star formers in the Universe were galaxies of stellar and 
gas mass of ~ IO^^Mq at redshifts 2 = 2 — 3, when the Universe was ~ 3 Gyr old. The 
common cases^'^ show star- formation rates (SFR) of 100 — 200 yr^^. These include 
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UV-selected galaxies termed BX/BM^^ and rest-frame optically selected galaxies termed 
sBzK/^ to be referred to collectively as "Star-Forming Galaxies" (SFGs). Their SFRs are 
much higher than the 4M0yr~^ in today's Milky Way, while their masses and dynamical 
times are comparable. The comoving space density of SFGs is n~2xl0~^ Mpc~^, implying 
within the standard cosmology (termed ACDM) that they reside in dark-matter haloes of 
masses ^3.5x lO^'^^Mo. The most extreme star formers are dusty Sub-Millimeter Galaxies 
(SMG)/^'i° with SFRs of up to ~ 1, 000 Mq jr'^ and n ~ 2 x 10"^ Mpc^l Whereas most 
SMGs could be starbursts induced by major mergers, the kinematics of the SFGs indicate 
extended, clumpy, thick rotating disks that are incompatible with the expected compact or 
highly perturbed kinematics of ongoing mergers. ^'^'^^'^ The puzzle is how massive galaxies 
form most of their stars so efficiently at early times through a process other than a major 
merger. A necessary condition is a steady, rapid gas supply into massive disks. 

It is first necessary to verify that the required rate of gas supply is compatible with the 
cosmological growth rate of dark matter haloes. The average growth rate of halo mass, M^, 
through mergers and smooth accretion, is derived^ based on the EPS theory of gravitational 
clusterings*^ (Supplementary Information, SI, §1), or from cosmological simulations. ^^'^^ For 
ACDM, the corresponding growth rate of the baryonic component is approximately 

M ~ 6.6 MiV' (1 + zf^' /.165 Mo yr-s , (1) 

where M12 = Mv/IO^^M©, and /.les is the baryonic fraction in the haloes in units of the 

cosmological value, /b = 0.165. Thus, at z = 2.2, the baryonic growth rate of haloes of 
2x 10^^ Mq is M ~ 2OOM0yr-\ sufficient for fueling the SFR in SFGs. However, the 
margin by which this is sufficient is not large, implying that (1) the incoming material 
must be mostly gaseous, (2) the cold gas must efficiently penetrate into the inner halo, and 
(3) the SFR must closely follow the gas supply rate. 

Penetrating Cold Ncirrow Streams 

The deep penetration is not a trivial matter, given that the halo masses of My > 
are above the threshold for virial shock heating,^^'^^'^'^ Mshock~ IO^^Mq. Such haloes are 
encompassed by a stable shock near their outer radius, R^, inside which gravity and thermal 
energy are in virial equilibrium. Gas falling in through the shock is expected to heat up to 
the virial temperature and stall in quasi-static equilibrium before it cools and descends into 
the inner galaxy. However, at 2; > 2, these hot haloes are penetrated by cold streams. ^'^ 
Dekel & Birnboim^ pointed out that because early haloes with > Mghock populate the 
massive tail of the distribution, they are fed by dark-matter filaments from the cosmic 
web that are narrow compared to and denser than the mean within the halo. The 
enhanced density of the gas along these filaments makes the fiows along them unstopable; 
in particular, they cool before they develop the pressure to support a shock, and thus avoid 
the shock heating (SI, §2). 

To investigate the penetration of cold streams, we study the way gas feeds massive high- 2; 
galaxies in the cosmological MareNostrum simulation — an adaptive-mesh hydrodynamical 
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Figure 1. Entropy, velocity and inward flux of cold streams pouring through hot haloes. The 
maps refer to a thin slice through one of our fiducial galaxies of = 1{)^'^Mq at z = 2.5. The 
arrows describe the velocity field, scaled such that the distance between the tails is 260kms~^. 
The circle marks the halo virial radius R^. The entropy, logif = log(T//9^/^), in units of the 
virial quantities, highlights (in red) the high-entropy medium filling the halo out to the virial 
shock outside R^. It exhibits (in blue) three, radial, low-entropy streams that penetrate into the 
inner disk, seen edge-on. The radial flux per solid angle is m = r'^pvr, in M0yr^^rad~^, where p 
is the gas density and Vr the radial velocity. 
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Figure 2. Streams in three dimensions. The map shows radial flux for the galaxy of Fig. 1 in a 
box of side length 320 kpc. The colours refer to inflow rate per solid angle of point-like tracers at 
the centers of cubic-grid cells. The dotted circle marks the halo virial radius. The appearance of 
three fairly radial streams seems to be generic in massive haloes at high redshift — a feature of 
the cosmic web that deserves an explanation. Two of the streams show gas clumps of mass on the 
order of one-tenth of the central galaxy, but most of the stream mass is smoother (SI, Fig. 10). 
The ^10^^ Mq clumps, which involve about one-third of the incoming mass, are also gas rich — 
in the current simulation only 30% of their baryons turn into starts before they merge with the 
central galaxy. 

simulation in a comoving box of side length 71 Mpc and a resolution of 1.4 kpc at the 
galaxy centers (SI §3). The gas maps in Figs. 1 & 2 demonstrate how the shock- heated, 
high-entropy, low-flux medium that fllls most of the halo is penetrated by three narrow, 
high-flux streams of low-entropy gas (SI, Figs. 7-10). The flux map demonstrates that more 
than 90% of the inflow is channeled through the streams (blue), at a rate that remains 
roughly the same at all radii. This rate is several times higher than the spherical average 
outside the virial sphere, mvir — 8M0yr~^rad~^ by eq. (4). The opening angle of a typical 
stream at -Rv is 20° — 30°, so the streams cover a total angular area of ~0.4rad^, namely a 
few percent of the sphere. When viewed from a given direction, the column density of cold 
gas below lO^K is above 10^°cm~^ for 25% of the area within the virial radius. While the 
pictures show the inner disk, the disk width is not resolved, so the associated phenomena 
such as shocks, star formation and feedback are treated in an approximate way only. 

The penetration is evaluated from the proflles of gas inflow rate, M(r), through shells 
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of radius r, Fig. 3 (SI, Fig. 11). The average profile reveals that the flow rate remains 
constant from well outside i?v~90kpc to the disk inside r~15kpc. 
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Figure 3. Accretion profiles M(r). Shown is the gas inflow rate through spherical shells of 
radius r, from the disk vicinity to almost twice the halo virial radius, obtained by integrating 
r'^pvr over the whole shell. The thick black curve is the average over the simulated galaxies of 
the fiducial case, ~ at z = 2.5. It shows deep penetration at a roughly constant 

rate ~lOOM0yr~^, consistent with the virial growth rate predicted by eq. (4). Apparently, the 
inflow rate does decay while traveling through the halo, but this decay is roughly compensated 
by the higher cosmological inflow rate when that gas entered the halo, eq. (4), leading to the 
apparent constancy of accretion rate with radius. The coloured curves refer to four representative 
galaxies, two showing clumps with ^ 0.1 and two with smoother flows involving only mini-minor 
clumps of /i < 0.1. Clumps with //^0.3 appear within 2i2v about once in ten galaxies; that is, 
major mergers are infrequent (SI, Fig. 11). The M(r) profiles serve for extracting the conditional 
probability distribution P(M|Mv), leading to the abundance n(> M) (SI, Fig. 12). 

Abundance of Gas Inflow Rates 

To relate the feeding by streams to the observed abundance of galaxies as a function of SFR, 
we use the MareNostrum inflow-rate profiles to evaluate n{> M), the comoving number 
density of galaxies with an inflow rate > M. We first extract the conditional probability 
distribution P{M\M^) by sampling the M(r) profiles uniformly in r, using the fact that 
the velocity along the streams is roughly constant (SI §§5,6). This is convolved with the 
halo mass function,^'' n{M^), to give 

n(M)= P{M\M,)n{M,)dM, . (2) 
Jo 
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The desired cumulative abundance n{>M), obtained by integration over the inflow rates 
from M to infinity, is shown at 2; = 2.2 in Fig. 4. Assuming that the SFR equals M, the curve 
referring to M lies safely above the observed values, marked by the symbols, indicating that 
the gas input rate is sufficient to explain the SFR. However, M and the SFR are allowed 
to differ only by a factor of ~2, confirming our suspicion that the SFR must closely follow 
the gas-input rate. Because at z^2.2 the star-forming galaxies constitute only a fraction 
of the observed ^10^^ Mq galaxies, ^^'^^ the requirement for a SFR SFR almost as great as 
M, based on Fig. 4, becomes even stronger. 
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Figure 4. Abundance of galaxies as a function of gas inflow rate, n{> M). Shown is the comoving 
number density, n, of galaxies with inflow rate higher than M at z = 2.2, as predicted from our 
analysis of the cosmological simulation. The upper curve refers to total inflow. It shows that 
galaxies with M > 150 Mq yr~^ are expected at a comoving number density n ~ 3 x 10^^ Mpc"'^ 
(similar to estimates in other simulations^ '''^^). Fluxes as high as M > 500 Mq yr~^ are anticipated 
at n ~ 6 X 10-5 Mpc"^. The lower curve is similar, but limited to gas input by /i>0.1 mergers. 
The symbols represent the vicinity of where the observed massive star-forming galaxies can be 
located once their observed SFR is identified with M. The sBzK/BX/BM galaxies are marked 
SFG.^^ The SMGs respectively brighter and fainter than 5 mJy are marked accordingly. ^^'^'^ We 
see that the overall gas inflow rate is sufficient for the observed SFR, but the small margin implies 
that the SFR must closely follow the rate of gas supply. Most of the massive star formers at a 
given SFR are expected to be observed while being fed by smooth flows rather than undergoing 
mergers. By studying the contribution of different halo masses to the abundance n(> M), we 
learn that the high-SFR SFGs and SMGs are associated with haloes of mass 10^^ — W^'^Mq (SI, 
Fig. 13). 
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Smooth Flows versus Mergers 

By analysing the dumpiness of the gas streams, using the sharp peaks of inflow in the 
M(r) profiles, we address the role of mergers versus smooth flows. We evaluate each clump 
mass by integrating Mdump = / M{r) dr/vr{r) across the peak, and estimate a mass ratio 
for the expected merger by /i = Mciump/(/b-^v), ignoring further mass loss in the clump on 
its way in and deviations of the galaxy baryon fraction from /b. We use the term "merger" 
to describe any major or minor merger of yU>0.1, as distinct from "smooth" flows, which 
include "mini-minor" mergers with fi<0.1. We find that about one-third of the mass is 
flowing in as mergers and the rest as smoother flows. However, the central galaxy is fed by 
a clump of /x>0.1 during less than 10% of the time; that is, the duty cycle for mergers is 
?7^0.1. A similar estimate is obtained using EPS merger rates^ and starburst durations of 
~50Myr at z — 2.5 from simulations^^ (SI, §5). 

Prom the difference between the two curves of Pig. 4 we learn that only a quarter of the 
galaxies with a given M are to be seen during a merger. The fact that the SPGs lie well 
above the merger curve even if the SFR is ~ M indicates that in most of them the star 
formation is driven by smooth streams. Thus, "SPG" could also stand for "Stream-Ped 
Galaxies". This may explain why these galaxies maintain an extended, thick disk while 
doubling their mass over a halo crossing time.^ On the other hand, if the SPR is ~M, we 
learn from Pig. 4 that about half of the bright SMGs and most of the fainter SMGs lie below 
the merger curve and are therefore consistent with being merger- induced starbursts.^^ 

Conclusion and Discussion 

We obtain that Stream-Ped Galaxies of baryonic mass ~10^^Mq at 2;~2.5 were the most 
productive star formers in the Universe. An integration of M over halo mass and time 
reveals that most of the stars in the Universe were formed in Stream-Ped Galaxies, within 
haloes >2 x lO^^M© at 1.5 < z < 4. The constraints on the overall SPR density at these 
epochs^° imply that SPR has been suppressed in smaller galaxies, e.g., by photoionization 
and stellar feedback. ^^'^^'^ The early presence of low-SPR galaxies^^'^^ requires quenching of 
SPR also at the massive end, perhaps due to gravitational heating by destructive streams. 

The streams are likely also to be responsible for compact spheroids, as an alternative to 
mergers^^ and the associated heating by expanding shocks. ^^'^^ Using eq. (4), we find that 
at z>2 the streams can maintain both the high gas fraction and the turbulence necessary 
for the disk to breakup into giant clumps by gravitational instability, with dispersion-to- 
rotation ratio a/V ~ 0.25, as observed.^ The clumps migrate inward and dissipatively 
merge into a spheroid. ^^'^^ The stream carrying the largest coherent fiux with an impact 
parameter of a few kiloparsecs determines the disk's spin and orientation, and the stream 
clumps perturb it. The incoming clumps and the growing spheroid can eventually stabilize 
the disk and suppress star formation. We can thus associate the streams with the main 
mode of galaxy and star formation occurring in massive haloes at 2; ~ 2 — 3; the streams 
that create the disks also make them fragment into giant clumps that serve both as the 
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sites of efficient star formation and the progenitors of the central spheroid, which in turn 
helps the streams to quench star formation. 

Although wet mergers may grow secondary disksSl, they are not as frequent as the 
observed SFGs (Fig. 4), these disks are neither gas rich nor clumpy enough, and, unlike 
most SFGs, they are dominated by stellar spheroids. 

The cold streams should be detectable by absorption or emission. For external back- 
ground sources, our simulation predicts that haloes with ~ IO^^Mq at z ~ 2.5 should 
contain gas at temperature < lO^K with cohimn densities > 10^°cm~^ covering ~ 25% of 
the area at radii between 20 and 100 kpc, with coherent velocities r$,200kms^^. Sources at 
the central galaxies should show absorption by the radial streams in ~ 5% of the galaxies, 
flowing in at ^200kms~^, with column densities ~10^^cm~^ (SI, Figs. 7-9). 
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SUPPLEMENTARY INFORMATION 

This is an extension of the Letter to Nature, aimed at providing further details, in support 
of the results reported in the main body of the Letter. 

1 Halo growth by EPS 

Neistein et al.^ used the EPS^^ theory of cosmological clustering into spherical haloes in 
virial equilibrium to derive a robust approximation for the average growth rate of halo 
virial mass My, 

dlnMv/du; = -(2/7r)V2[a2(Mv/g)-(72(Mv)]-^/^ uj = S,/D{t). (3) 

The time variable cu, which makes the expression time invariant, is inversely proportional to 
D{t), the linear growth rate of density fluctuations at time t in the assumed cosmology, with 
Sc ^1.68. The power spectrum of initial density fluctuations enters via the variance a'^{My). 
The constant q is 2.2 with an uncertainty of ±0.1 intrinsic to the EPS theory. Eq. (3) has 
been confirmed to resemble the assembly rate in cosmological A^-body simulations.^^ 

For the ACDM cosmology,"^^ a flat Universe with 72% dark energy, mass dominated 
by cold dark matter, and fluctuation normalization parameter erg = 0.8, the corresponding 
growth rate of the baryonic component is well fitted by the practical formula^ 

M ~ 6.6 M^y^ (1 + zf-^' /,i65 Mq yr-i , (4) 

where = M^/IO^^Mq, and /. 165 is the baryonic fraction in the matter assembled into 
haloes in units of the cosmological value /b = 0.165. 

2 On the origin of nairrow streams 

Dekel & Birnboim 2006^ (hereafter DB06), following the simulations by Keres et al. (2005, 
§6.2) and their own analysis of the simulations by A. Kravtsov, pointed out that at redshifts 
higher than Zcrit ~ 2, narrow cold streams penetrate deep into the dark-matter haloes even 
when the haloes are more massive than the shock- heating scale, Mghock ~ 10^^ Mq, and 
proposed a possible explanation for this phenomenon. This prediction is summarized in 
Fig. 5. 

The critical condition for a stable virial shock is that the radiative cooling rate behind 
the shock is slower than the compression rate, ^cooi <^comp) allowing the buildup of pressure 
support behind the shock against global gravitational collapse. Based on a spherical anal- 
ysis, DB06 found that a virial shock should exist in dark-matter haloes above a threshold 
mass Mghock ~ IO^^Mq that is rather constant in time, at an actual value that is sensitive 
to the metallicity of the gas in the halo. The existence of such a threshold mass and its 
value as a function of redshift have been confirmed by the analysis of cosmological simu- 
lations. ^'^'^^'^ However, at high redshifts, even above the threshold mass, a shock is not 
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Figure 5. Analytic prediction for the regimes dominated by cold flows and shock-heated 
medium in the plane of halo mass and redshift, based on Fig. 7 of DB06. The nearly 
horizontal curve marks the robust threshold mass for a stable shock based on spherical 
infall analysis, Mshock(^)- Below this curve the flows are predicted to be predominantly 
cold and above it a shock-heated medium is expected to extend out to the halo virial 
radius. The inclined solid curve is the conjectured upper limit for cold streams, valid at 
redshifts higher than Zcrit ~ 2. The hot medium in haloes of > Mghock at z > Zcrit is 
predicted to host penetrating cold streams, while haloes of a similar mass at z ^ ^cv\t 
expected to be all hot, shutting off most of the gas supply to the inner galaxy. Also shown 
is the characteristic Press-Schechter halo mass M^.(z); it is much smaller than Mghock at 
z>2. 
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expected to develop along narrow, cold, radial streams that penetrate through the halo, 
because the cooling there is more efficient than in the surrounding halo. 

The appearance of intense streams at high z, as opposed to their absence at low z, is 
likely to reflect the interplay between the shock-heating scale and the independent charac- 
teristic scale of nonlinear clustering, i.e., the Press-Schechter^^ mass M^, that corresponds 
to the typical dark-matter haloes forming at a given epoch. The key difference between the 
two epochs is that the rapid growth of with time, as seen in Fig. 5, makes Mshock3>M* 
aX z>2 while Mshock ~ M^, at lower redshifts. 

Cosmological A^-body simulations^'^^ reveal that while the rare dark-matter haloes of 
My ^ tend to form at the nodes of intersection of a few filaments of the cosmic web, 
the typical haloes of My ~ tend to reside inside such filaments. Since the filament width 

1 /3 

is comparable to the typical halo size R^ocM^ and seems not to vary much with position 
along the filament, one expects the rare haloes to be fed by a few streams that are narrow 
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compared to the halo size, while the typical haloes accrete from a wide angle in a practically 
spherical pattern. Assuming that at any given epoch the accretion rate of dark matter, M, 
is roughly proportional to the halo mass (Eq. 1 of the Letter), while the virial densities 
in haloes of all masses are the same (by definition), the geometrical difference implies that 
the densities in the filaments penetrating 3> haloes are higher by a factor of a few 
than the typical densities in their host haloes. The above is demonstrated in Fig. 6 (as well 
as in Fig. 5 of Ocvirk, Pichon & Teyssier 2008). 




Figure 6. The pattern of dark-matter inflow in a shell (1— 3)-Rv outside two sample haloes 
from a cosmological N-body simulation at 2; = (based on P. Seleson & A. Dekel, in 
preparation). Left: a typical halo with ~ M*. Right: a rare halo with ^ M^. 
In terms of the different ways by which these two haloes are fed by dark-matter, they 
correspond to two haloes of the same mass ~ lO^^M©, but at z ~ and z ~ 2 — 3 
respectively. Top: dark-matter density contrast about the mean density in the shell. 
Middle: particle velocity dispersion (loosely termed "temperature"), in terms of the virial 
value. Bottom: infall velocity, represented by the anisotropy parameter /3, where (3 = 
corresponds to isotropic velocities and /3 = 1 to pure radial motions. We see that the typical 
halo resides inside a broad filament so it is practically fed by wide-angle diffuse accretion. 
On the other hand, the rare halo is fed by narrow, dense, radially in-fiowing filaments. 

Assuming that the density of the gas flowing along the filaments scales with the dark- 
matter density, and that the infall velocity is comparable to the halo virial velocity, one 
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expects the cooling rate in the filaments feeding an My 3> halo to be higher by a factor 
of a few than in the surrounding spherical halo. If the compression rate in the filaments 
is comparable to that in the host halo, this implies that the thin filaments should have 

a harder time supporting a stable shock. As a result, the critical halo mass for shock 
heating in the filaments feeding it must be larger by a factor of a few. This is the case for 
Mv ^ Mshock haloes at high redshifts. 

A crude estimate along the lines above led DB06 to the conjectured upper limit for 
penetrating streams shown in Fig. 5: 

Mstream ~ Mghock for /M* < Mghock , (5) 

where the characteristic width of the streams is oc (/M*)^/"', with / a factor of order a few. 
At low z, where /M* > Mghock, cold fiows exist only for < Mshock- At high z, where 
/M* < Mshock, cold streams appear even in Mv>Mshock haloes where shocks have heated 
part of the gas, as long as My < Mstream- The critical redshift ^^crit separating these two 
regimes is defined by 

/M*(Zcrit) = Mshock . (6) 

This crude maximum mass for cold streams is shown in Fig. 5 for an arbitrary choice of 
/ = 3. 

A preliminary analysis of the MareNostrum simulation^ confirms the crude prediction 
of eq. (5), when taking into account the lower metallicities in the simulation compared to 
that assumed in the analytic calculation. The streams analysed in the current Letter, in 
dark-matter haloes of My — IO^^Mq at 2; = 2.5, represent an encouraging confirmation 
of the validity of the DB06 conjecture. Further analysis in progress (T. Goerdt et al., in 
preparation) indicates, for example, that at z = 2.5, the fraction of infiow in cold streams 
drops by a factor of three at My ^ 2 x IO^^Mq, much in the spirit of the crude prediction 
of Fig. 5. The permitted cold gas supply by streams in massive haloes at high redshift, 
followed by a shutdown above Mshock at low redshifts, turn out to provide good match 
to many observed galaxy properties when these features are incorporated in semi-analytic 
simulations of galaxy formation. ^^"^^ Still, the dependence of the stream properties on 
redshift and halo mass is yet to be explored in a more quantitative way. 



3 The MareNostrum simulation 

The cosmological simulation^'"^^ used in the present analysis has been performed with the 
Eulerian AMR code RAMSES'^'^ on 2,048 processors of the MareNostrum supercomputer. 
The code simulates the coupled gas and dark-matter dynamics, using a Particle-Mesh 
scheme for the dark-matter component and a second-order Godunov scheme for the gas 
component. In order to describe the formation of dense star-forming disks, the code includes 
metal-dependent radiative cooling, UV heating by a standard photo-ionizing background, 
star formation, supernovae feedback and metal enrichment. The simulation box of co- 
moving 50/i~^Mpc involved 1,024^ dark- matter particles and 4 x 10^ gas cells. Using a 
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quasi-Lagrangian refinement strategy, the spatial resolution reaches ~ 1 h kpc in physical 
units at all times. The dark- matter particle mass is ~ 10''' Mq, so each of the haloes studied 
here consists of ~ 10^ particles within the virial radius. Since one can rehably describe the 
formation of haloes down to ~ 100 particles, namely ~ IO^Mq, the IO^^Mq haloes ad- 
dressed here are three orders of magnitude above the minimum halo mass. This simulation 
allows us to capture the important properties of gas accretion into galaxies in more than 
100 haloes of ~ 10^^ M© at 2;~2.5, thus providing a large statistical sample. A first analysis 
of galaxies from this simulation^ have confirmed the bi-modal nature of cold fiows and hot 
media as a function of mass and redshift.^ 

Our current analysis is based on robust features that are properly simulated, such as 
the large-scale structure of the streams, the total flux in them, and the gas clumps larger 
than IO^Mq. However, the finite resolution does not allow a fair treatment of small- 
scale gas phenomena such as turbulence in the hot gas, ram-pressure stripping of clumps, 
hydrodynamical instabilities at the stream boundaries, and the formation of small clumps. 
Furthermore, the current resolution does not allow a detailed study of the disks that form 
at the halo centers as the disk thickness is barely resolved. More accurate analysis of the 
fine stream structure and disk buildup should await simulations of higher resolution. 

4 Maps of entropy, flux and density 

Figs. 7 to 9 extend the visual information provided by Figs. 1 and 2 of the Letter. They 
display different gas properties that highlight the structure and kinematics of the cold 
streams in three simulated galaxies of M^ — IO^'^Mq at z — 2.5. 

The entropy maps show log(T/p^/^) where the temperature and gas density are in units 
of the virial temperature and mean density within the halo virial radius R^. They exhibit 
the virial shock, covering most of the area of the virial sphere and sometimes extending 
beyond 2R^. The narrow streams are of much lower entropy, by more than three orders of 
magnitude, comparable to the low entropy in the central disk they lead to. The boundaries 
between the streams and the hot medium within the virial radius are sharp and well defined. 
We also note that semi-cylindrical shocks sometimes partly surround the elongated streams 
long before they enter the halo virial radius. 

The arrows mark the velocity field projected on the slice plane, and the fiux colour 
maps show the flow rate per solid angle, m = r^pvr- The flux inward is almost exclusively 
channeled through the narrow streams, typically involving 95% of the total inflow rate. 
This flux is several times the average over a sphere, rhyir — 8 yr~^rad~^. The opening 
angle of a typical stream at Ry is 20 — 30°, so the streams cover a total area of ~0.4rad^, 
namely a few percent of the sphere. The velocity fleld in the hot medium is turbulent and 
sometimes showing vast outflows. The inward flux over most of the sphere area is negligible, 
both inside and outside the virial radius or the virial shock. The streaming velocities are 
supersonic, with a Mach number of order a few. 

Although the streams tend to be rather radial when viewed on scales comparable to 
the halo virial radius, some of them flow in with impact parameters on the order of 10 kpc. 
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comparable to the disk sizes. The steady high flux along a line of a rather fixed orientation 
with a non-negligible impact parameter is the source of angular momentum required for 
the buildup of an extended rotating disk (A. Zinger et al., in preparation). 

The gas density maps emphasize the narrowness of the streams, and reveal that they 
are typically denser than the surrounding medium by more than an order of magnitude. 
This confirms the prediction described in §2, and explains why a virial shock is avoided 
along the streams, allowing them to penetrate cold and unperturbed into the inner halo. 

The column-density maps of the in-flowing material are obtained by summing up the 
densities in grid cells along each line of sight inside the box of side 320 kpc. The cells 
that enter this sum are only those where the inward flux per solid angle is at least twice 
the average over a sphere based on Eq. 1 of the Letter. These maps highlight the three- 
dimensional configuration of radial streams, and the clumps along some of them. Such 
column density maps will serve us in producing detailed predictions to guide observations 
in search of the cold streams in high-redshift galaxies, in absorption and in emission. 

Fig. 10 displays three-dimensional TIPSY* pictures of the radial influx m, similar to 
Fig. 2 of the Letter. It shows the overall structure of the in-flowing streams in 3D perspec- 
tive for four simulated haloes. The pictures reveal that the typical configuration is of three 
major narrow streams. Some of the streams are straight lines, and others are curved. Some 
of the streams are of rather fixed width from well outside R^, and others display a conical 
shape, starting broad at large radii and getting narrower as they penetrate into the halo. 
The gas streams show dense clumps, with about one third of the stream mass in clumps 
of mass ratio > 0.1, namely mass above ^lO^^M®. The rest is in smaller clumps, some 
clearly hidden below the resolution limit. Since these mini-minor clumps are not expected 
to cause significant damage to the central disk,^^ we can refer to them in this respect as 
"smooth" flows. It is not clear at this point to what extent the smooth component is truly 
smooth or built by mini-minor clumps, and whether the perfect smoothness has a physical 
origin or is merely a numerical artifact, but this distinction does not make a qualitative 
difference to our present diskussion. 



* http : / /www- hpcc. astro .washington.edu /tools /tipsy /tipsy.html 
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Figure 7. Gas in lialo 314 of the MareNostrum simulation. Three maps refer to a thin 
equatorial slice. They show (a) entropy logK = log(T/p^/^) in units of the virial quantities, 
(b) radial flux m = r'^pVr in yr~^ rad~^, and (c) log density in atoms per cm~^. The 
fourth panel shows log column density through the 3D box of side 320 kpc, in cm~^, con- 
sidering only the cells where the radial flux inward is at least twice as high as the average 
over a shell based on Eq. 1 of the Letter. The circle marks the virial radius. 
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Figure 8. Gas in lialo 303 of tlie MareNostrum simulation. See Fig. 7. 
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Figure 9. Gas in lialo 311 of tlie MareNostrum simulation. See Fig. 7. 
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Figure 10. Inward flux in the three-dimensional boxes of side 320 kpc centered on galaxies 
314, 303, 311 and 310 from the MareNostrum simulations. The colours refer to inflow rate 
per solid angle of point-like tracers at the centers of cubic-grid cells. The dotted circle 
marks the virial radius. All haloes show high-flux streams, some smooth and some with 
embedded clumps. Galaxy 310 (bottom right) is undergoing multiple minor mergers due 
to the particularly clumpy streams. 
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5 Accretion profiles and probabifity distribution 




R [kpc] 

Figure 11. Profiles of total gas inflow rate through spherical shells as in Fig. 3 of the 
Letter. Shown here are twelve galaxies of ~ IO^^Mq at z = 2.5, randomly chosen from 
the simulation. The lower six panels show clumps that correspond to mergers of mass ratio 
fi > 0.1, while the upper six are fed by smoother flows with only mini- minor mergers of 
fi < 0.1. 

Fig. 11 is an extension of Fig. 3 of the Letter, presenting the influx profiles of twelve 
galaxies, all with ^ IO^^Mq at z = 2.5, chosen at random from the MareNostrum 
simulation. The profiles extend from r = 15 kpc, the disk vicinity, to r = 160 kpc, almost 
twice the virial radius of -Rv — 90kpc. 

In order to evaluate the conditional probability P{M\My) that enters Eq. 2 of the Letter, 
we first measure Po{M\Mq) from a fair sample of MareNostrum haloes of our fiducial case 
Mq = lO^^M© at zq = 2.5. This probability distribution, shown in Fig. 12, is derived 
by sampling the M(r) profiles shown in Fig. 11 uniformly in r, using the fact that the 
inflow velocity along the streams is roughly constant. The tail of the distribution shown 
in Fig. 12, at M > 200MQyr~^, is dominated by yU > 0.1 mergers, while the main body 
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Figure 12. The conditional probability distribution P{M\My) for the fiducial case 
M^ = Mo = W^Mq at z = 2.5. 

of the distribution is mostly due to smoother streams. Recall that the average is about 
lOOMgyr-^ 

We then generalize Po(-^l-^o) to other masses using the scaling from eq. (4), M oc 
M^'^^, namely 

P(M|M.) = Po[M(Mo/M.)^-^^|Mo] . (7) 

At 2;~2.5, this scaling of M is good to within a factor of two for < IO^^Mq, beyond 
which it becomes a more severe overestimate (Goerdt et al., in preparation). The results 
for other redshifts {z>2) are obtained using the scaling from eq. (4), M oc (1 + z)"^'"^^. 

Half the galaxies shown in Fig. 11 turn out to show clumps leading to mergers of /i> 0.1, 
and the rest show only smaller clumps in smoother fiows. One can read from the relative 
width of the clumps in Fig. 11 that the duty cycle for /i> 0.1 clumps in each individual 
galaxy is less than 0.1. By comparing the areas above the individual profiles with the 
average for galaxies of that mass and redshift, one can see that on average only about one 
third of the stream mass is in clumps. 

A similar estimate is obtained from EPS, by reading from Fig. 6 of Neistein &; DekeF the 
rate dN/duj of >yU mergers into a halo M^. The typical starburst duration is At~0.1-Rv/K' 
(~50 Myr at 2; = 2.5), based on merger simulations^^ or the typical peak width in the M(r) 
profiles (Fig. 11), given streaming at the virial velocity = {GM^/RyY^'^ ~220kms~^. 
This leads to ri= {dN/ dcu) At 0.09 for M = 2xW^Mq at z = 2.2. 
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Figure 13. Comoving number density of galaxies witli total gas inflow rate higher than 
M at 2; = 2.2, as in Fig. 4 of the main body of the Letter. The numbers in ma- 
genta next to the right axis refer to log of haloes with the corresponding abundance. 
Left: Dependence on the maximum halo mass that contributes to cold streams, for 
M,toeam = 10^^ 10^2.5^ ^0^3, IO^^Mq. Right: Variation with redshift, z = 2,3,4,5,7. 



6 The abundance as a function of mass and redshift 

As described above, the conditional probability distribution P{M\M^) has been estimated 
by scaling the results from the simulated haloes of = lO^^M©. Preliminary analysis 
of more massive haloes at that redshift (T. Goerdt et al., in preparation) indicates that 
the actual inflow rate starts dropping below the adopted estimate in haloes more massive 
than Mstream ~ IO^^Mq. For a first crude estimate of the effect this might have on our 
results shown in Fig. 4 of the Letter, we re-compute the comoving number density n{>M) 
as described in the main text, but now limit the halo mass range that contributes to M 
by an upper cutoff at Mstream- Fig. 13 shows the results for different values of Mstream- 
We see that a cutoff at Mstream = IO^^Mq makes only a small difference to n{> M), by a 
factor of ~ 2 at the high-M regime corresponding to the bright SMGs. Thus, the decay 
of cold streams above IO^^Mq is not expected to alter our results in a qualitative way. On 
the other hand, we learn from the fact that the symbols for SFGs and bright SMGs lie 
far above the lower curve that the high-SFR objects at these redshifts are dominated by 
central galaxies in haloes more massive than IO^^Mq. In fact, we read from the figure that 
some of the SFGs and many of the bright SMGs are associated with haloes more massive 
than 3 x 10^^ Mq. 

Fig. 13 also shows the predicted abundance n(>M) at different redshifts, now applying 
no finite upper mass cutoff Mstream- This is justified for z > 2 based on our preliminary 
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investigation of the MareNostrum galaxies at different redsliifts and masses and consistent 
witfi tfie conjecture of DB06^ shown in Fig. 5. We see that the comoving abundance of 
galaxies with M~ 150 Mq yr~^ is predicted to vary by a factor less than two between z — 2 
and 4. By 2; ~ 7 that abundance drops by an order of magnitude. The variation with redshift 
is somewhat larger at the high-fiux end, toward M~10^MQyr~^. At lower redshifts, the 
contribution of streams in massive haloes above Mshock is most likely overestimated by 
this procedure, so a similar analysis in the low- 2; regime should impose an upper limit at 

-^^stream — -^^shock ■ 



